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Abstract 

Soft limits of inflationary correlation functions are both observationally relevant and theoretically 
robust. Various theorems can be proven about them that are insensitive to detailed model- 
building assumptions. In this paper, we re-derive several of these theorems in a universal way. 
Our method makes manifest why soft limits are such an interesting probe of the spectrum of 
additional light fields during inflation. Wc illustrate these abstract results with a detailed case 
study of the soft limits of quasi-single-field inflation. 
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1 Introduction 



Soft limits of primordial correlation functions are of special phenomenological relevance, since 
many observational probes of non-Gaussianity — such as scale-dependent halo bias [1] and CMB 
//-distortions [2] — involve a separation of scales. At the same time, soft limits are theoretically 
clean. Focusing on soft limits allows one to prove a number of interesting theorems that are 
insensitive to detailed model-building assumptions, but are sensitive to the number of light degrees 
of freedom during inflation. Moreover, as we will explain in this paper, the momentum scaling 
in the soft limits is sensitive to details of the spectrum of particles during inflation. 

There are two classes of soft limits: 

1. The 'squeezed limit' of an A^-point function refers to taking one external momentum to 
be smaller than all others (see Fig. 1). Maldacena's consistency relation for single- field 



inflation [3] (see also [4-6] ) relates the squeezed limit of the three-point function of curvature 
perturbations C to the variation of the two-point function under dilation. Similar results 
have subsequently been derived for A-point functions [7-10]. These results are explained 
by realizing that the curvature perturbation ( non-linearly realizes a dilation symmetry in 
any FRW background [11] (see also [12-14]). Any violation of the consistency relations 
would rule out all single-field models and point to additional degrees of freedom during 
inflation. 

2. The 'collapsed limit' of an A-point function takes one internal momentum (i.e. the sum 
of M of the A external momenta) to be smaller than all external momenta (see Fig. 2). 
The Suyama-Yamaguchi inequality [15] (see also [16, 17, 20]) puts a bound on the collapsed 



limit of the four-point function in terms of the squeezed limit of the three-point function. 
At tree level, the inequality is saturated if a 'single source' creates the primordial curvature 




Figure 1: Soft external momentum or the 'squeezed limit'. 
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Figure 2: Soft internal momentum or the 'collapsed limit'. 
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perturbation and its non-Gaussianity [21]. Collapsed limits are therefore probes of 'multiple 
source' (or 'hidden sector') non-Gaussianities. An interesting observational signature of 
this effect is stochasticity in the scale-dependent halo bias [22] (see [23] for another possible 
signature) . 

A boost in the collapsed limits arises naturally in models of quasi-single field inflation 
(QSFI) [24]. In these models a light inflaton field mixes with an additional scalar field (the 
'isocurvaton') with mass of order the Hubble scale H} This mixing allows large non-Gaussianities 
in the isocurvaton sector to be communicated to the observable sector. The time dependence 
of the massive field on superhorizon scales leads to a characteristic scaling behavior of A^-point 
functions in their soft limits [25]. Moreover, the perturbative coupling between the observable 
sector and the isocurvature sector leads to the boost in the collapsed limit relative to the squeezed 
limit. 

The outline of the paper is as follows: In Section 2, we present new proofs of various 
important theorems [3, 15] concerning the soft limits of primordial correlation functions. Although 
these results aren't new, some of our methods and perspectives may be. In particular, we make 
a bit more precise the analogy between soft limits of inflationary correlation functions and soft 
pion limits of QCD [26-31]. In Appendix A, we present a number of technical details and 
generalizations related to the theorems of Section 2. In Section 3, we compute the soft limits of 
QSFI explicitly in terms of the fundamental parameters of the theory. We conclude in Section 4 
with a discussion of the prospects of using non-Gaussianity as a particle detector. 



2 Generalities: Soft Limits of Inflationary Correlation Functions 

We begin with proofs of various important theorems [3, 15] concerning the soft limits of infla- 
tionary correlation functions. 

2.1 Soft Internal Momenta 

We first consider the collapsed limit of the four-point function, i.e. we let the internal momentum 
of the four-point function become soft, while keeping all external momenta hard. We will show 
that the four-point function in this limit is bounded from below by the squeezed limit of the 
corresponding three-point function. Our proof will make the role of additional degrees of freedom 
manifest. It will also make obvious how the proof generalizes to higher A'^-point functions. 

2.1.1 Suyama-Yamaguchi Bound 

Definitions. — We characterize the soft limits of the three- and four-point functions by the follow- 
ing two quantities 

^ ^ 1 li^ (Sh^k^ , (2.1) 
•' 12fci^o P1P2 

= 4 t^™n PPP ' (^-^^ 

^Such a spectrum of particle masses arises generically in supersymmetric theories of inflation [25]. 
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where (•)' denotes the expectation value without the factor {2-k)^6{Y^ ki). Moreover, we have 
defined Pi = P({ki) = {C^ Cj^.)' and ki2 = \ki + ^2!- All correlation functions are to be understood 
as 'connected' correlation functions. For local non-Gaussianity, the definitions (2.1) and (2.2) re- 
duce to the familiar (momentum-independent) definitions of /nl [32] and t^l [21]. More generally, 
however, /nl and Tnl can be momentum dependent — e.g. in §3, we will find /nl(^i) oc kf^"' and 
t'nl(^12) oc A;^^^", with 3 > q > |. We will now give a proof of the bound t^l > (|/nl)^ [15, 17]. 
In fact, we will give a proof of a slightly more general inequality, which includes Tnl > (f /nl)^ 
as a special case. 

Collapsed four-point function. — We start by considering the following four-point function,^ 

{cHx,u)e{o,t.)) = m\x,t,)cHo,um) . (2.3) 

where is some time at which (, has become frozen and is the vacuum'^ of the interacting 
theory. Next, we insert a complete set of momentum eigenstates \n^) in the correlation function, 

(C'(x,t.)C'(0,t.)) = ^(C'(f,t.)|n,-)(ng.|C'(0,t.)) , (2.4) 

n,q 

where Ylnqi') = Sn/g(") ~ Yin I (27rp (")' Working in the Heisenberg picture and 

therefore the states are formally^ independent of time. Here, q is the total momentum of a 
state and n is a label for individual states (e.g. 'particle number'^ states or different particle 
species) within the momentum eigenspace. In the following, we will drop the time label t^,. Using 
(C^(x)|ng') = e*'^'^(C^(0)|nq-), we can write eq. (2.4) as 

(C'(x)C'(O)),^ = / d'xe-''-^C\x)C\6)) = |(n,HC'(0))|2 . (2.5) 

Notice that the r.h.s. is written as a sum over positive definite terms. If we can show that one of 
the contributions is the square of the squeezed limit of the three-point function of then we can 
conclude that the collapsed four-point function — the l.h.s. of (2.5) — is always greater or equal. 

Squeezed three-point function. — To complete the proof, we insert a complete set of states in the 
three-point function {C{x)C'^{0)) , i.e. 

= J d3xe+*-^(C(f)C'(0)) = ^ (C,Hn^>(«dC'(0)) , (2.6) 

n,q 

where = |Cg)^^ = (f^|C_^- It is always possible to choose the basis states \n^), such that 

{C^\n^) = {2^f5{k - q) 6in Pl'^ik) . (2.7) 

— 1/2 

Up to an arbitrary phase, this implies |1^) = P^ (k) iCj:)- In an abuse of the standard particle 
physics terminology, we will call |lg) the 'one-zeta' state. All other states, which we will denote 

^Our conclusions will be unchanged if we replace ("^(x, f^) — >■ (^{x, ti,)(^{x + e, t*) and C^(0, t*) ^ C(0, t*)C(e*, i*), 
which may be useful in order to avoid divergences (or contact terms) when the operators are at coincident points, 
''it is worth remarking that our results hold even if is not the vacuum state, as long as the state is invariant 
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vacuum 
« ICfe) "one-zeta" 
\nj:) _L ICj) "everything else" 

Table 1: Spectrum of states and their meanings. 

by {nj:), are orthogonal to They may be associated with 'single-particle' states corresponding 
to additional degrees of freedom or simply 'multi-zeta' states, i.e. states created with higher powers 
of the operator (. Substituting (2.7) into (2.6), we get 

(gc'(o)) = (i,HC'(o))p^'/'(A:). 

Bound on the collapsed four-point function. — We can therefore write eq. (2.5) as 



(2i 



PAk) 



+ Ei^^fcic'(o))p 



(2.9) 



Since K"'fclC^(0))P ^ 0, we conclude that 



(C'(x)c^(o)), > 



KCglC^(o))l^ 



(2.10) 



If we take the soft limit A; — t- 0, we get a relation between the collapsed four-point function and the 
squeezed three-point function. Moreover, from our derivation it is easy to see how the result can 
be generalized to give a relation between the soft limits of 2A^-point functions and the squeezed 
limits of (A^ + l)-point functions, 

(C^(f)C^(o)),^ = ^^C^ + ElKHC^(o))l^ . (2.11) 

Suyama-Yamaguchi bound. — Eq. (2.10) can be written as 



> 



qi-k' 



91 92 



which, in the limit /c — )• 0, becomes 



f^^Pc{qi)Pdq2) > 
qi q2 qi 52 



(2.12) 



(2.13) 



under spatial translations and rotations. 

''in practice, we define the Hilbert space perturbatively using the free theory. The states then become implicitly 
time dependent through the projection onto the interacting vacuum. 

^Of course, the notion of 'particle states' is subtle in time-dependent cosmological backgrounds without asymp- 
totically free states. None of those subtleties will affect our arguments, so we will, for convenience, sometimes 
(ab)use the particle physics terminology. 
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So far this is completely general. To get the specific form of the Suyama-Yamaguchi bound [15] one 
has to make a further assumption. It is often the case that and (/nl)^ are either momentum 
independent or have the same momentum dependence (so that we can still extract a momentum- 
independent amplitude). We then get 



rNL>(i/NL)' . (2.14) 



2.1.2 Interpretation 

What is the physical meaning of the extra contributions in eq. (2.9)? It is well-known [15, 21] 
that curvature fluctuations arising from a 'single source' (which may or may not be the inflaton) 
saturate the Suyama-Yamaguchi inequality (at least at tree level). In this case, the squeezed 
three-point function of C completely determines the collapsed four-point function — i.e. the extra 
contributions in (2.9) are subdominant, in the sense that 

lim , ^ , = (single source at tree level) . (2.15) 
*^^o (IglC^) 

Non-negligible (j^-^^j^qIC^) can have the following origins: 

1. Extra fields. A straightforward way to have contributions from extra states is to have 
multiple sources for the primordial fluctuations. This can lead to a significant boost of the 
collapsed four-point function relative to the square of the squeezed three-point function, 
"Tnl ^ (|/nl)^- As a measure of such 'multi-source' non-Gaussianity we can define 



X = lim 



(C^(f)C^(o))gPc(fe) ^ 



limVW))^. (,,6) 



We will see an explicit example of this situation in Section 3. 

2. Gravitons. Strictly speaking, even single-source inflation always has additional states cor- 
responding to the graviton, i.e. tensor metric perturbations. In fact, in single-field slow-roll 
inflation these states give the dominant contribution to the collapsed four-point function, 
Tnl ~ 0{e) [10], where e = -H/H'^ < 1. This is parametrically larger than f^^ ~ C(e^) [3], 
although both contributions are too small to be observable. Being interested in observable 
signatures, we will typically suppress the graviton contributions in our results. It would be 
easy to include their effects as 0(e) corrections. 

3. Excited initial states. The effects of excited initial states on primordial non-Gaussianity have 
been explored in [35-40]. In particular, it was pointed out [37] that non-trivial contributions 
to soft limits are possible. It is conceivable that these effects are captured by (?^^^q|C^) ii^ 
our formalism, but we leave a detailed analysis for the future [41]. 

4. Loops. Even in single-source inflation we can (formally) have t^l > (|/nl)^ from loops [18- 
20]. In single-clock inflation [33], these loop contributions aren't physical and can be re- 
moved by the right choice of physical coordinates [34] . In more general single-source scenar- 
ios, loops may contribute non-trivially. Having said that, in this paper, we won't consider 
these loop effects, but instead focus on tree- level contributions. 
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2.2 Soft External Momenta 

We now use similar ideas to derive the famous consistency relations for correlation functions with 
soft external momenta [3, 4]. Our approach will draw inspiration from methods used in proving 
the low-energy theorems of theories with spontaneous symmetry breaking [29, 30], such as the 
classic soft pion theorems of QCD [26-28]. This connection was emphasized recently in [11]. The 
proofs in this section will be restricted to the limit /c — )• 0, but we will show in Appendix A how 
they can be generalized to finite k. 

2.2.1 Symmetries of Inflation 

The consistency relations arise from the symmetries of inflation, which we therefore quickly review 
(for further discussion see [11, 42, 43]). 



4>{t) 




Figure 3: Conformal diagram of de Sitter space. The inflationary baci<ground picks out a preferred time-slicing. 

Inflation is characterized by a quasi-de Sitter background 

ds2~-di2 + e2^*dx2 . (2.17) 

The isometry group of de Sitter space, SO{A, 1), contains 6 translations and rotations, 1 dilation 

Ht ^ Ht- X , X ^ {1 + X)x , (2.18) 

and 3 special conformal transformations (SCTs) 

Ht ^ Ht-2b-x , X ^ f +6(F"^e~2^*-x^) + 2(6-f)x , (2.19) 

where A and b are infinitesimal parameters. At late times, t — t- oo, these isometrics act as 
conformal transformations on the spatial boundary 

X ^ (1 + A)x, (2.20) 
X ^ X -hx'^ + 2{h- x)x . (2.21) 

The time-dependence of the inflationary background (j){t) spontaneously breaks the symmetry 
under dilations and special conformal transformations. Inflationary correlations are therefore. 
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in general, not invariant under the conformal transformations. However, since the symmetry is 
non-hnearly reahzed, we still expect non-trivial relations between different A^-point functions. 

As explained nicely in recent work by Creminelli, Noreha, and Simonovic [12] and Hinter- 
bichler, Hui and Khoury [11], the symmetries in eqs. (2.20) and (2.21), in fact, have significance 
beyond quasi-de Sitter backgrounds. They are symmetries of adiabatic modes in any FRW back- 
ground.^ The Goldstone boson associated with the spontaneous breaking of the symmetry is the 
curvature perturbation^ C, which in co-moving gauge is defined by [47, 48] 

mjix,t) = a\t)e^'^^^''^5ij . (2.23) 

As we will see, the squeezed limit of A^-point functions of C are related to the variation of (A^ — 1)- 
point functions under dilations and special conformation transformations. This is a consequence 
of the conformal transformations being non-linearly realized. 

We now show how these considerations can be incorporated into our formalism. For con- 
creteness, we focus on the dilation symmetry^ in co-moving coordinates 

X ^ {1 + X)x , (2.24) 
C ^ C + A(l + x • (9^C) • (2.25) 

This symmetry is a large^*^ diffeomorphism that is non-linearly realized by the "gauge-invariant" 
perturbation ^. It should not be confused with an isometry of FRW (which is linearly realized) 
or a residual^^ gauge transformation. During inflation, only time diffeomorphisms are broken 
by the background [33]. Therefore, this large spatial diffeomorphism will act linearly on all the 
non-gravitational degrees of freedom. Associated with the dilation symmetry of the action is a 
conserved current and a charge Qd = J d^x J^{x,t). Although the charge does not exist after 
symmetry breaking (due to an IR divergence), local commutators are still well defined, e.g. 

6dC = i[Qd,C] = -l-x-dsC ■ (2.26) 

For future reference, we will need the matrix element {il.\QdCq\^) = {Qd\Cq)- Invariance of 
the vacuum under spatial translations and rotations, fixes the form of this matrix element to be 

{Qd\Q) = if{q) ■ i27TfS{q) , (2.27) 

^The symmetries are still explicit for fields that don't couple to the inflaton, such as tensor metric fluctua- 
tions [42] or scalar spectator flelds [43, 44]. Even for inflaton fluctuations, the symmetries are restored in the 
so-called decoupling limit, where Mpi — > oo, — > 0, with M^iH — > const. 

^Adiabatic modes are deflned locally as coordinate redeflnitions of the unperturbed FRW solution. Coordinate 
transformations of the form (2.20) and (2.21) with time-dependent parameters \{t) and b{t) induce the following 
co-moving curvature perturbation 

Q = \{t) + 2b{t) ■ X . (2.22) 

For constant transformation parameters A and b this gauge mode is the fc — > limit of a physical solution [11, 12], 
a la Weinberg [45] , with suitable fall-off at spatial inflnity. 

*For spontaneously broken 'spacetime' symmetries the number of Goldstone bosons does not have to equal the 
number of broken generators [46] . 

^It is our understanding that a complete treatment will appear in [13, 14] (see also [12]). 

^"A "large diffeomorphism" is one that does not vanish at infinity. These transformations are distinct from the 
gauge redundancies of the theory, namely local diffeomorphisms that vanish at infinity. 

^^There is a residual gauge transformation in "^-gauge", eq. (2.23), that does not act on [49]. 



8 



where / is a function that depends only on the magnitude of the 3- momentum q = \q\- In order 
for this to be consistent with the commutator (2.26), we require Re[/(g)] = ^. The imaginary 
part of f{q) won't contribute in what we are about to say, so we can set it to zero. We then have 

(Qdig = h ■ C^^f^iQ) ■ (2-28) 

2.2.2 Consistency Relations 

We now use these results to give a quick proof of the 'Adler zero' [26] for the squeezed limit of 
the three-point function in single-field inflation [3, 4]. We will also see how extra fields have to 
enter in order to produce a non-trivial signal in the squeezed limit. 

Associated with the conserved current is the following Ward identity [50] 

i {J^{x, t)C{y, QC{z, )) = d{t - U)S{x - y) {ddCiy, U) C{z, U)) 

+ 5{t - U)5{x - z) (C(y, U) ddCiz, U)) , (2.29) 

where for clarity we have this time shown all the time dependences explicitly. If we integrate^^ 
eq. (2.29) over dt / d^x, we arrive at 

([Qd, C^^C^J) = ([Q., CuM) + ^Cfei t^-^' ^ajJ) ' (2-30) 

where we have reverted to dropping the common time coordinate t^. In going from (2.29) to 
(2.30), we have dropped the integral over diJ\ because it only receives contributions from spatial 
infinity (see Appendix A). With the help of (2.26), we can express the r.h.s. of (2.30) as 

{[Qd. CkSk,\) = -^(3 + ■ d^^^ + h ■ d^^JiC-^^CkJ' {27TfS{ki + h) . (2.31) 

Here and in the following, we only keep terms that contribute to the connected part of the 
correlation functions. Next, we insert a complete set of states in the l.h.s. of (2.30), 

iiQ<i^C,;SkJ) = E[^^'^l^^)KHCfe^.,)-(QA^J",-)(^d^3^)] > (2-32) 

n,q 

where 

{Qd\ng) = ^Cn{q) ■ {27rfS{q) . (2.33) 

Since cq = {Qd)' = {^\Qd\^)' £ the vacuum insertion doesn't contribute. We separate the 
'one-zeta' state, with ci{q) = ^^'^{q) (up to a removable phase), from the rest^'^ 



i2Trf 6{h + k2) . (2.34) 



^^The integral over time is extended in the imaginary direction because our correlation functions are "time 
ordered" with respect to the integration contour in the complex plane (see Appendix A for further details). 

^''The functions |cn(q)|^ in eq. (2.34) are analogous to the spectral density functions p(g^) in the Kallen- 
Lehmann representation of correlation functions [29]. In principle, one could determine the properties of |cfi(g)p 
from {J^{x)j;^id)). 
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In deriving eq. (2.34) we have used the fact that 

where the second equality comes from the invariance of correlation functions under parity. Com- 
paring eqs. (2.31) and (2.34), we get 



(2.36) 



This result generalizes straightforwardly to the squeezed limits of higher A^-point functions (see 
Fig. 1) 



(CglCfc^ • • • Cfc^) g-s>0 



N 



3(iV-l) + 5]^..9gJ(Cg^...Q^)' 

^Im[cfi(g)(n,-|Cfc, •••Cfc^)' 



i=l 



(2.37) 



Moreover, it is easy to see how this result implies consistency relations for soft internal mo- 
menta [51] (see Fig. 2), 



PdQ) 



Pdi) Pdi) 

L E(Cgi • • • C4j-.-)'("dC.^,,,, ■ • • C,^,)' > (2-38) 



Pd<l) 



where q = ^f^i ki. Again, this was derived by inserting a complete set of states and separating 

— 1/2 

the 'one-zeto' state 11^-) = {'l)\Cq) from all states orthogonal to it, \ng-). 

2.2.3 Interpretation 

There are two important consequences of the q ^ limit. First, in this limit causality explains 
the absence of a potential boundary term arising from the diJ^ term in (2.29) (see Appendix A 
for the complete argument). Second, in this limit the 'multi-particle' states among \h^) are not 
expected to contribute the leading divergent term, so that we can drop them on the r.h.s. of (2.36). 
It is straightforward to check this in most examples, however, a general proof is beyond the 
scope of this paper. 

In the special case of single-field inflation, eq. (2.36) therefore reduces to the famous single- 
field consistency relation [3] 

(CglCfc^Cfcj) g->0 



PdQ) 



3 + ^1 • % + ^2 • dj:^ 



iCkSk,y = {i-ns)Pdki) . (2.39) 



'^^Even in the context of S-matrix elements for the scattering of Goldstone bosons, there are known exceptions 
to this statement [29] (§19.2). 

^^Excited initial states may contribute corrections at finite q, or when the number of e-folds is taken to infinity 
before taking g — >■ (see e.g. [36-40]). This is related to the observation that q ^ fci_2 is not a sufficient condition 
for satisfying the consistency relation (see e.g. [6, 54, 55]). Understanding these finite q effects can be important 
given that observations will only probe a finite range of momenta. 
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For scale-invariant fluctuations/^ i.e. (CgjC^g)' ~ (^1^2)"^''^ or Ug = 1, the r.h.s. of (2.39) is zero 
and the squeezed hmit of the three-point function vanishes. This is analogous to the Adler zero 
of soft pion physics [26] (see also [29, 30]). 

In multi- field inflation, we can get contributions to the squeezed limit, even for scale- invariant 
fluctuations, 

^ -El">[=.-.(«)<«.HCE.Cfe>'] . (2.40) 

This expression doesn't just tell us that we need extra flelds to generate a squeezed limit, but also 
how the extra fields must enter. We will study an explicit example of this situation in Section 3. 

3 Case Study: Soft Limits of Quasi-Single-Field Inflation 

In this section, we will illustrate some of the abstract results of the previous section with the 
concrete example of quasi-single-fleld inflation (QSFI) [24]. This is a particularly interesting case 
study, since QSFI exhibits non-trivial behaviour in its soft limits: 

i) the soft limits of A^-point functions have characteristic momentum scalings that depend on 
the mass of the isocurvaton field; 

ii) the collapsed limits of 2A^-point functions are naturally boosted. 

In §3.1, we explain these features of QSFI qualitatively. In §3.2, we compute the soft limits of 
the three- and four-point functions explicitly. 

3.1 Review of Quasi-Single-Field Inflation 



Figure 4: Particle spectrum of quasi-single-field inflation. 



Dynamics. — QSFI couples a massive 'isocurvaton' field a to an 'adiabatic' perturbation 9 (see 
Fig. 4). The dynamics of these fields is governed at leading order by the following free field 
Hamiltonian (density) 

no = l{d^9f + ^{d^af + lm''a\ (3.1) 

^^Scale-invariance is natural if the theory has an internal shift symmetry, (l> (f> + c. Both this symmetry and 
the 'spacetime' dilation symmetry of de Sitter space are spontaneously broken by the time-dependent inflationary 
background, (j}{t). However, a linear combination of shift and spacetime dilation is preserved [52]. At late times, 
this diagonal symmetry becomes the spatial dilation symmetry that we discuss in this section. This explains the 
scale- invariance of the two-point function [.53]. Softly breaking the shift symmetry allows for weakly scale-dependent 
fluctuation spectra. 
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where m ^ H. At linear order, the adiabatic perturbation and the curvature perturbation <^ 
are related by 

e = -V2eM^iC, (3.2) 

where e = —H/H^. A quadratic mixing term, e.g. ^mix = —pGcr, converts a into 9 (or C,). This 
can communicate large interactions in the isocurvaton sector to the observable sector. Our basic 
example will be the cubic interaction ^3 = fia^ . Following [24], as part of the 

interaction Hamiltonian 

'Hint = -pOa + /xfj^ . 
In order to describe the mixing perturbatively we require p < H. 
Quantization. — The fields are quantized in the interaction picture 



(3.3) 



%(r) = Ufc(r)ag + 4(r)a[- , 

where dr = dt/a{t) is conformal time and lifc(r) and Vkir) are the free-field mode functions 

H 



Uk{r) 



2k^ 
n H 



(1 + ikT)e 



-ikr 



, (-fcT)3/2H(;)(-A;r) , with v = \r- -"^ . 



(3.4) 
(3.5) 

(3.6) 
(3.7) 



Here, H[/^' is the Hankel function of the first kind. The massive mode decays on superhorizon 
scales, t'fc(T) ~ (— r)^/^^'^. This feature completely explains the non-trivial scaling of correlation 
functions in the soft limits [25]. The mode functions have been normalized such that [dg, at] = 

[6^, ftt ] = (27r)^(5(A; -|- A:'), which makes d^ (dt) and 6g (bt) the standard annihilation (creation) 

operators. Since there is no d-a coupling in T^Qj we also have [dg, St] = 0. We compute correlation 
functions of products of field operators, Q{t), in the in-in formalism [55, 56] 



(Q(r)) = {n\ 



fexpii dr'Fint(r') 



Qir) 



Texp 



dr'i?int(r') 



\n) , (3. 



where the symbols T and T denote time-ordering and anti-time-ordering, respectively. The 
interaction Hamiltonian Hi^t = f d^x a'^ T-Lmt can be written as a sum of two terms 



^2 = - j d^xa^pO'a and = j 



(3.9) 



The mixing term H2 is treated as a perturbative insertion in Feynman diagrams (see Fig. 5). 
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Figure 5: Perturbative description of Af-point functions in quasi-single-field inflation. 

Soft limits. — We end this review section with a brief summary of the soft hmits of QSFI: 
• The three-point function has the following momentum scaling in the squeezed limit 



lim M^M. oc kt" . (3.10) 

This result has a simple explanation [25]: The long (or soft) u-mode exits the horizon at 
the time ri ~ k^^ , after which it decays with a rate that depends on its mass, akir) ~ 
(— r)^/^"*^. Its correlation with the short (or hard) modes is evaluated at the horizon crossing 
of the short modes, i.e. at T2 ~ /c^^. By this time, the amplitude of the long cr-mode has 
decayed by a factor of {t2/ti)'^^'^^'^ = (^1/^2)^/^"''. This explains the suppression of the 
squeezed limit in eq. (3.10). Below we will confirm this intuition with an exact calculation. 

The four-point function in the collapsed limit scales as 

hm — V, ri \ K2 ■ (3-11) 

fci2^o Pciku) 

Again, the scaling follows from the superhorizon evolution of the soft ci-mode. 



~ » ilUf ■ (3.12) 



The amplitudes of the collapsed four-point function and the squeezed three-point function 
are related as 

We can understand this from the Feynman diagrams in Fig. 5. While in the three-point 
function every cr-propagator attaches to an external 9 (or Q, in the four-point function the 
virtual cr-mode does not have to be converted to C. The four-point function therefore only 
gets four powers of the small coupling p/H, while the square of the three-point function 
has six powers of p/H. For fixed three-point amplitude, we therefore get a boost in the 
four-point amplitude by a factor of {p/H)~'^ » 1. 

3.2 Explicit 'in-in' Computations 

We now want to go beyond the qualitative discussion of the previous section, and compute the 
soft limits of the three- and four-point functions explicitly in terms of the fundamental parameters 
of the theory. 
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3.2.1 Power Spectrum 

To assess the regime of validity of perturbation theory, we first use the in-in master formula (3.8) 
to compute the leading correction to the power spectrum of the adiabatic mode 9, 



\ fci fc2 ' 

Substituting eqs. (3.6), (3.7) and (3.9), we find [24] 

1 H"^ " 



k^PAk) 



pi I I 



i + .M(|) + 



where 



c{v) = 27rRe 



dxi 



drc2 ( x^'/'hW (xi)e^"ix-'/'H(,^' (x2)e 



-1X2 



1^-1/2tt(2), 



10 
8 
fi 



(3.13) 



(3.14) 



(3.15) 



J . . , . . , — ^ p 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 

The function c{v) formally diverges in the limit i/ — )• |. In this limit, the isocurvaton a becomes 
massless and freezes after horizon exit, instead of decaying. As a result, the transfer of cr to C 
continues for an indefinite time after horizon exit and evaluating the curvature perturbation at 
the infinite future r = is no longer justified. Instead, the bispectrum should be evaluated at 
the time Tf when inflation ends. This regulates the divergence. Finally, we remark that the 
calculation is under perturbative control if and only if 

c(^) ■ (3.16) 

3.2.2 Soft Limits: Methodology 

We now derive the soft limits of QSFI in a style that will be familiar to cosmologists [4]. At the 
end of §3.2.3 and §3.2.4, we will demonstrate the precise equivalence between this method and 
the formalism of §2. 

Long-short split. — In cosmology, soft limits are described as the correlation between long and 
short modes. Following [3, 4, 58], we treat the long- wavelength fields as c/assim/ backgrounds for 
the short-wavelength quantum fields. The fields 9 and a are split into long and short modes, 

9 = 9i^ + 9s and cr = + as, (3.17) 

^^For the specific mixing term we consider, T^mix ~ —p6a, an effective mass of order p is generated for a [24, 57]. 
If we were to re-sum tlie perturbation series, we expect tliat tlie divergence to be regulated by tlie effective mass. 
However, it isn't guaranteed tliat this possibility exists for generic mixing interactions, so we will ignore it here as 
well. 

^^This is morally similar to the split into soft and hard modes in the factorization theorems of QCD [59]. 
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where 



^3 k 



(T)e 



ik.x 



and 



k>A 



(2^)3 



T e 



Ak.x 



(3.18) 



ik<\ (27r)' 

and similarly for a. The precise value of the momentum cut-off A won't be important. We then 
write the interaction Hamiltonian H\^t in terms of 0^, 6s-, Cl and as- The quadratic mixing 
term H2 doesn't produce any coupling between long and short modes. The only relevance of the 
mixing is therefore to convert (7$ into 9s via 



-pa' 



j d^x Use's 



(3.19) 



The cubic interaction H^, on the other hand, does couple long modes to short modes. The only 
interaction that will be relevant in the soft limit is 



Hr 



3/ia'* / d'^XCTifj 



(3.20) 



Soft limits. — Our goal is to use eqs. (3.19) and (3.20) to compute the squeezed limit of the 
three-point function, 



and the collapsed limit of the four-point function, 

(^ki^k2^kjkj ^^^^ ((%i%2)'^i(%3^fc4^'^i^ • 

We see that both soft limits can be expressed in terms of the modulated power spectrum 



(3.21) 
(3.22) 

'ki^'ki'^'L > (3.23) 

i.e. the power spectrum {^^^Oj^J in the presence of a long- wavelength background (Tj,. We compute 
eq. (3.23) in the in-in formalism 



H-,^k.2l''L{r) 



dri 



dT2 



dT3( Fint(T3), HiAT2), ^int (n ) , %^ %^ (r 



where H\^t = Hj^ss + -f^ss- In the limit r — t- 0, this becomes 



kik2''^L 



{A12 + B12) 



(3.24) 



(3.25) 



2 m kf 

where the ^-term is obtained by replacing -ffint(''"3) by Hj^ss and the two other interaction terms 
by Hss, while the i?-term is obtained by replacing -ffint(7"2) by Hj^ss and the rest by Hss- This 
leads to 



A12 = 
B12 = 



dxi 



—00 




XI 



dX2 



dxi 



—00 

Xl 



X2 



dX2 



— 00 
X2 



dx3 fA{xi,X2,X3) {crL)k^j^k^{x^/ki) 
dxa /b(xi,X2,X3) (o-i,)g^_^^^(x2/A;i) 



(3.26) 
(3.27) 



with 

fA{xi 



{-xiX2{x:if)-^'^Bhi{-xi)sin{-X2) X ImfH(f)(-a;i)H(J^n-^2)(H«(-X3))' 



fB{xi) = (-xi(x2)'x3)-^/2sin(-xi) X Im(H(^'(-xi)H(^'(-X2))lm(HW(-X2)H(^'(-X3)e^"'3) . 
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3.2.3 Bispectrum in the Squeezed Limit 

With this preliminary work, it is now straightforward to compute the squeezed hmit of the 
bispectrum, 



fci \ k2 k'^i'^Ll 



2 fci 



A 



23/ 



All 



(3.28) 



The only slight complication is due to the fact that the calculation requires the unequal time 
correlation function of two fields, i.e. (fT^(T)^g'(O)). However, since 0^{t) freezes in the limit 
Igrl <C 1, we have lim O^It) OffU)) and hence 



(^,^(r)%(0)) {'J^,{r)e^{r)) , 



(3.29) 



where 



(^fc(r)W) = -i I df / a(f)V(0|<Tp-(f)e'_p-.(f)a,KT)^,-.(r)|0) + c.c. 

'OO J p 



= -2p{27rf 6 (k + q) Re 

Using the super horizon limits of the mode functions, 

H 



i / dfa{ffvk{f)vl{T)u {f)u*{T) 



lim Ukir) 



and lim Vkir) 

|fcr|«l 



we get 



where 



(cT,(r)0,-.(r)) 



|fcr|«i /2 r(z^) p 



IT 2"^-^ H fc3 



-/fcr)3/2-'^ d(i/) (27r)35(p + g) 
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d{v) 



/ dx(-2;)-^/2j^e [H«(-x)e" 

J —oo 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Substituting (3.32) into (3.28), we obtain the squeezed limit of the bispectrum for the curvature 
perturbation C,, 



ikSkM.)' ^ f/NL(-)-Pc(^i)^c(fc2)(^ 



k2) 



3/2-y 



where 



/,,(z.)^p(i.)d(i.).^.A-i 



(3.34) 
(3.35) 
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Here, d(z/) is the same function as in eq. (3.33) and p{u) is defined as 

rO 



T(i>) r' /"^'^ 

= 5^3/2 \_l ^^1 / 



{-X^f'^-'fAiXi) + {-X2f'^''' fB{x^) . (3.36) 



Eq. (3.34) precisely matches the result that one obtains by taking the limit fei — )• of the full 
expression of the bispectrum [24] . 

In the following insert we show how to reproduce eq. (3.34) from our analysis in §2: 



We start by defining a complete set of one-particle states 



\lf)^P^"\i) and \lf)^P;'f^\aij 



In order for the states to be orthogonal, we have defined 
As our expansion parameter we define 



< 1 . 



(3.37) 



(3.38) 



(3.39) 



To determine the effect of a on the squeezed limit of the three-point function of C, we first compute the 
parameter cg- in eq. (2.40). Since spatial dilations act linearly on a, we have 



This implies that 



caiq) 



.(Cfk?)' 



PrPi 



1/2 



(3.40) 



(3.41) 



The coefficient cg. may have a non-zero real part, but it won't contribute in what follows. We therefore 
get 



Im 



c.(9)(iriC'^' 



PcPa 



Im [z(a,-|C')'] 



PcPa ^ 



and eq. (2.40) becomes 



Pdl) 



Pd<i)PAq) 



(3.42) 



(3.43) 



This is equivalent to the starting point of this section, eq. (3.28), so following the same steps will lead to 
the answer in eq. (3.34). 
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3.2.4 Trispectrum in the Collapsed Limit 

The collapsed limit of the four-point function is equal to the correlation of two modulated power 
spectra, 



2 J kfkl 



(3.44) 



where the sum over contractions is 

w = (n/'^'^^) (n/ '^^^) ^^-^^^ 

+ /B(x.)/A(y.)(K),^^^(f )K),^3^(||)> + /i.(^.)/i.(y.)(K),^,,(t)K)..3,(t)) 

Using 

22-lr(zy)2 



Kl(n)fT,-2(r2)) 
this becomes 



91,92- 



r(^) 



{-qirif/^~''(,-q2T2f^-' (27r)3<5 (gl + q^) , (3.46) 



22-21-^ (A;iA;3)3/2-^ (A;i2)2^' 



(2^)3 5(A?i 



ny"dx,^ (/^(3;,)(-X3)=^/2-.^^^(^^)(_^^)3/2-. 



Eq. (3.44) can therefore be written as 

A:i2-s-0, 



/i4,,2 174 1 

H'^J- ^6 (feife3)9/2- (A:i2)2- 



(3.47) 



(3.48) 



where 



500 
400 
300 
200 
100 



(3.49) 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 



p 



^ The divergence for 1/ — t- arises from the 



The function t{v) diverges as — t- and v , 2 
fact that the expansion of the Hankel function in eq. (3.31) isn't justified when ki2/k <C e~^/'^ 
Instead the expansion of the Hankel function should be changed to 



(3.50) 
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This regulates the divergence. The divergence of t{u) for — )• | has the same origin (and the 
same resolution) as the divergence of c(z^) in (3.15). Finally, converting (3.48) to curvature 
perturbations, we get 



3-21/ 



where 



TNL(i^) 



(|/nl(^))^ 

e2(^) 



with e{i') = — d{i>) . 

H 



(3.51) 



(3.52) 



Note that we have recovered the scaling in eq. (3.11). Since (i^(z^) ~ c(i^), the condition on 
perturbative control, eq. (3.16), becomes 



(3.53) 



Hence, as long as the calculation is under perturbative control, we find Tnl > (|/nl)^) consistent 
with the Suyama-Yamaguchi inequality [15]. 

Finally, we show how eq. (3.51) also follows from our analysis in §2: 

We insert the complete set of one-particle states, eq. (3.37), into the four-point function 



(c)= J [(ciif)(i^^'io+(ciir>(irio 



where 



11^"' 



N 



\l'^:>)-a\l'^> 



with a = (1^'|1^')' e{iy) and N = 1 + 0{s). At leading order in e, eq. (3.54) becomes 



(3.54) 



(3.55) 



(3.56) 



This is equivalent to the starting point of this section, eq. (3.44), so following the same steps will lead to 
the answer in eq. (3.51). 
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4 Discussion: Non-Gaussianity as a Particle Detector 



single-field inflation 



pi 



H - 



quasi-single-field inflation 



multi-field inflation 

A 



lim k^lQCt Cr ) 














» 




> 


J/nl) 



Figure 6: Particle spectra of inflationary models and their observational signatures. 

Observations of the CMB and LSS probe the physics of inflation through the statistics of the 
primordial curvature perturbations. Soft limits of the correlation functions of these perturbations 
are especially sensitive to the number of light degrees of freedom during inflation. In this paper, 
we have explained how the spectrum of particles can impact observations (see Fig. 6). 

Single-field inflation. — The fluctuations in single-field inflation are Goldstone bosons of a spon- 
taneously broken dilation symmetry [11]. Being non-linearly realized, the broken symmetry is 
still reflected in Ward identities relating the squeezed limits of A^-point functions to the sym- 
metry variation of (A'^ — l)-point functions [11, 42, 43]. Up to small corrections associated with 
the breaking of scale-invar iance, the squeezed limits of all A^-point functions therefore vanish in 
single-fleld inflation [3]. This is analogous to the Adler zero [26] of chiral perturbation theory [29]. 
In this note, we have given a new proof of this famous result. Measuring a significant three-point 
function in the squeezed limit would therefore automatically rule out all single-field models. Go- 
ing beyond single-field infiation, the squeezed limit captures information about the spectrum of 
additional light degrees of freedom. 

Multi-field inflation. — Coupling additional fields to the inflaton relaxes some of the symmetry 
constraints on inflationary correlation functions. The extra fields can therefore produce a sig- 
nificant signal in the squeezed limit. In fact, if the non-Gaussianity arises from superhorizon 
evolution, then the signal tends to be peaked in the squeezed configuration. If a single field 
(but not necessarily the inflaton) creates the primordial curvature perturbations ( and their non- 
Gaussianity, then the collapsed limit of the four-point function is uniquely fixed by the squeezed 
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limit of the three-point function. However, this correlation is broken if more than one field sources 
C and its interactions. In that case one can show that the collapsed four-point function is always 
larger than the square of the squeezed three-point function [15]. In this note, we have given a 
new proof of this result. 

Quasi-single-field inflation. — Supersymmetric models of inflation often mix a light inflaton field 
with an additional isocurvaton field whose mass is of order the inflationary Hubble scale [25]. Such 
models of quasi-single-field inflation [24] naturally give rise to a boosted four-point function and to 
characteristic momentum scalings in the soft limits. We have computed these signatures explicitly 
in terms of the fundamental parameters of the theory. As we have explained elsewhere [22, 25], 
these results are relevant for connecting the microscopic physics of QSFI to late-time observables 
such as stochastic halo bias [22] . 

The above shows that soft limits of primordial non-Gaussianity provide a remarkable op- 
portunity to probe the spectrum of particles during inflation, i.e. at energies totally inaccessible 
to conventional particle physics experiments. The PLANCK satellite [60] will measure the CMB 
anisotropics over a large range of scales and hence will significantly improve the observational con- 
straints on the soft limits of primordial non-Gaussianity. At the same time, the Sloan Digital Sky 
Survey [61] continues to collect data on the clustering of collapsed objects. The scale-dependent 
biasing [1] of these objects is a very sensitive probe of the squeezed limit of the primordial corre- 
lation functions [62]. Looking ahead, LSS-experiments like EUCLID [63] and LSST [64], as well 
as CMB-experiments like PIXIE [65] and CMBPol [66] all have the potential to teach us a great 
deal about the soft limits of inflationary correlations. 
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A Consistency Relations at Finite Momentum 



The soft limits described in §2.2 were presented in the strict g — )• limit. However, since these 
results follow from current conservation alone, one might expect that they can be extended to 
finite q. In this appendix, we will demonstrate that this is indeed the case. 

By definition, any conserved charge is written as an integral of the current over a spatial 
slice. Throughout the main text, we have worked directly with the charge Q = J dc^x J^{x,t) to 
derive our results. In particular, we have considered the following commutator 

iiiQXimm = {6QC{y)c{z)) + m)SQC{z)) , (a.i) 

where it is understood that all fields for which the time dependence isn't shown are evaluated at 
the same time t^,. Eq. (A.I) follows from current conservation via the Ward identity [30, 50] 

i ( J^(x, t)C{y, U)C{z, U,))t = Sc{t- U)5{x - y) {5^^ U) ({z, U)) 

+ 6c{t - U)6{x - z) {C{y,U)6QC{z,Q) , (A.2) 

where the terms proportional to 6cit — t^,) arise from the time derivative acting on the "time- 
ordering" operator T. This feature of the Ward identity deserves a few more words of explanation: 
By definition, the above in-in correlation functions are "time-ordered" on a contour C that runs 
from t = —oo + ie to — oo — ie, passing through t^. Formally, we can write this as 

{Oi{ti)02{t2))T = 0{Xi - As) (Ol(tl)02(t2)) + e{X2 - Ai) {02{t2)Oi{h)) , (A.3) 

where 6 is the Heaviside function and A is a real number that parametrizes the contour C — e.g. 
we may choose X = t — t^, on the "incoming" contour C+ and A = — t on the "outgoing" 
contour C_ (see Fig. 7). The delta function 5c{t — t*) on the r.h.s. of (A.2) is therefore defined as 

A f ^ / '^(*-**) onC+, 

Oc[t-ti.) = < ... ^ (A. 4 

I —d[t^ — t) on C_ . 

If we now integrate (A.2) over the contour jl*_^l_^l^dt J d^x = <f dt J d^x through t = t^, we 
arrive at eq. (A.I). We have dropped the integral over 9j J* because it only receives contributions 
from spatial infinity. 




t 



C- 

Figure 7: Illustration of the contour employed in the computation of in-in correlation functions. 

When translated into momentum space, the above treatment is limited to g = 0. However, 
nothing required that we integrate the Ward identity in this way. Alternatively, we can multiply 
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eq. (A. 2) by e^'^'^ and then integrate over x and t, 

dt I d'x e^-?- i5(-)(J^(f,t)C(y)C(^"'))T = e^'^-y{6QCmm) + e''^-^({y)5QC{z)) . (A.5) 
Taking the hmit e, At — )• and going to momentum space, we find 

Inserting a complete set of states in the first term, one arrives at the identity 



n,p 



(A.7) 



Note that this result holds without making any assumptions about the relative or absolute sizes 
of q, ki and k2- This isn't surprising, since so far we have only used current conservation which 
holds at all momenta. In order to arrive back at eq. (A.l), we must show that we can indeed 
neglect the extra term in eq. (A.7) in the limit g — t- 0. The choice of contour in the complex time 
domain implies that 

dt q,{Jl{t)Cj;^Cf:,) ^ ^dtq,{[Jt{t),Cj:^C^:.J) . (A.8) 

Causality requires that [J*(x, t*), Civ^ **)C(-2*5 i*)] vanishes when x is outside the light-cone of y 
and z. For x ^ y,z, the operators are space-like separated and the commutator is a proportional 
to a delta function in position space, or in momentum space [</^,Cfc] = /('i')^'^^^-, where f{q) is 
an analytic function. Translation invariance ensures that the correlation function (O^^^^^j:^) can 

be written in terms of /c2. Therefore, the final term in (A.8) will be analytic in q. Furthermore, 
invariance under spatial rotations requires that the commutator vanishes at g = (using the 
rotation ki —ki). These two facts together imply 

lim {[Ji{t), C^^CkJ) = fih, k2) q' . (A.9) 

Hence, we find 

lim ^ [{4\n^) (n,HC,^^C,^,) - h.c] = (5^C,-.+,^, C,^,) + (C.^ 6^C,^^^,J + 0{q') . (A.IO) 

n,p 

The conclusion that the consistency relation holds at finite q up to corrections suppressed by q^ 
is extremely important. It ensures that contributions to the l.h.s. of (A.IO) that vanish like 
for a > 1 cannot be mimicked by finite q effects in single-field infiation. See [6] for an alternative 
derivation of this result. 



The causal structure of the background will be modified by metric fluctuations. In all situations of interest in 
this paper, the fluctuations are small and this effect is negligible [34]. On the other hand, in eternal inflation the 
presence of large metric fluctuations can modify the causal structure signiflcantly [67]. 
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